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Abstract 

The Brownie tag-recovery model is useful for estimating harvest rates but 
assumes all tagged individuals survive to the first hunting season; otherwise, 
mortality between time of tagging and the hunting season will cause the Brownie 
estimator to be negatively biased. Alternatively, fitting animals with radio trans- 
mitters can be used to accurately estimate harvest rate but may be more costly. 
We developed a joint model to estimate harvest and annual survival rates that 
combines known-fate data from animals fitted with transmitters to estimate the 
probability of surviving the period from capture to the first hunting season, and 
data from reward-tagged animals in a Brownie tag-recovery model. We evaluated 
bias and precision of the joint estimator, and how to optimally allocate effort 
between animals fitted with radio transmitters and inexpensive ear tags or leg 
bands. Tagging-to-harvest survival rates from >20 individuals with radio trans- 
mitters combined with 50-100 reward tags resulted in an unbiased and precise 
estimator of harvest rates. In addition, the joint model can test whether transmit- 
ters affect an individual's probability of being harvested. We illustrate application 
of the model using data from wild turkey, Meleagris gallapavo, to estimate harvest 
rates, and data from white-tailed deer, Odocoileus virginianus, to evaluate whether 
the presence of a visible radio transmitter is related to the probability of a deer 
being harvested. The joint known-fate tag-recovery model eliminates the require- 
ment to capture and mark animals immediately prior to the hunting season to 
obtain accurate and precise estimates of harvest rate. In addition, the joint model 
can assess whether marking animals with radio transmitters affects the individ- 
ual's probability of being harvested, caused by hunter selectivity or changes in a 
marked animal's behavior. 



Introduction 

Individually marking wildlife is a common technique used 
to obtain data for estimating population parameters, such 
as survival and harvest rates. In some cases, more than 
one data type is available for an individual animal, such 
as recaptures, resightings, telemetry locations, or recover- 
ies. Recent advances in mark-recapture models include 
joint models, in which a single model incorporates infor- 
mation from multiple data types. Such models can 
increase the number of estimable parameters, provide 
greater precision of estimates, and reduce estimator bias 



(Nasution et al. 2001). Models have been developed for 
joint analysis of many combinations of data types (Burn- 
ham 1993; Lebreton et al. 1995; Barker 1997; Catchpole 
et al. 1998; Powell et al. 2000; Nasution et al. 2001; Bar- 
ker et al. 2004; Gopalaswamy et al. 2012). A combined 
radio-telemetry and tag-recovery model was developed 
(Pollock et al. 2004) and implemented (Bacheler et al. 
2009) to estimate instantaneous harvest and natural mor- 
tality rates within a fisheries context. Although Lakhani 
(1987) used radio-telemetry data to inform age-specific 
annual survival rates from a program of marking only 
juveniles of a nongame species, a model that incorporates 
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radio-telemetry into a tag-recovery estimator for terres- 
trial game species has not been developed. 

Radio -telemetry and tag-recovery programs have strengths 
and weaknesses associated with implementation and data 
collection that have to be evaluated in light of research or 
management goals. Monitoring individuals via radio telem- 
etry can provide more detailed information about individ- 
ual fates, and these data can be used to estimate survival 
with a Kaplan-Meier known-fate estimator (Kaplan and 
Meier 1958). For example, year round data on survival, 
cause of death, dispersal, and habitat relationships can be 
collected for specific individuals (Kenward 2001). Con- 
versely, it can be costly to monitor a sufficient number of 
individuals fitted with radio transmitters to detect biologi- 
cally important differences among groups (Powell et al. 
2000; Kenward 2001; Winterstein et al. 2001; Rogers and 
White 2007). Also, individual fates must be monitored for 
the duration of the study and cost depends on the number 
of field technicians, size of the study area, number of 
marked individuals, and how frequently an individual must 
be located (Kenward 2001). 

In addition, hunters may notice the presence of a radio 
transmitter before deciding to harvest an individual 
(Mayer et al. 2002; Jacques et al. 2011); however, it is 
unclear if detecting the transmitter influences their deci- 
sion to attempt harvest (Fuller 1990; Etter et al. 2002; Jac- 
ques et al. 2011). Some hunters may not harvest 
individuals with a radio transmitter because they believe 
that it is illegal (Jacques et al. 2011) or they are benefit- 
ting the research program by foregoing harvest (F. E. 
Buderman, pers. obs.). A reluctance to harvest an animal 
with a radio transmitter would result in underestimating 
the harvest rate of the population. However, if hunters 
select animals with radio transmitters because of the nov- 
elty, this would result in overestimating the harvest rate. 
In addition, there may be concern that the transmitter 
itself leads to behavior by the marked individual that 
makes it differentially susceptible to harvest (CasweU et al. 
2012). 

In comparison with a radio-telemetry study, a study 
capturing and marking individuals with inexpensive tags 
(hereafter called "tags") should have lower costs per ani- 
mal marked, although additional expenses are incurred 
when reward tags from marked animals are subsequently 
reported. Consequently, more animals can be marked in a 
tagging study for the same total cost as a radio-telemetry 
study, although information is only available on harvested 
individuals (Brownie et al. 1985). A problem with the 
Brownie tag-recovery estimator, however, is the assump- 
tion of no mortality between the time of tagging and the 
first recovery period (e.g., hunting season). For example, 
white-tailed deer, Odocoileus virginianus, are most easily 
captured in winter, when there are few alternative food 



sources and individuals congregate (e.g., Verme 1973; 
Whitlaw et al. 1998; Norton et al. 2012), but not har- 
vested until autumn; mortality during the 6-9 months 
between tagging and harvest, though minimal, may nega- 
tively bias the harvest rate estimator. Diefenbach et al. 
(2012) successfully used a Brownie model to estimate 
spring harvest rates of male wild turkeys, Meleagris gallop- 
avo, but such a model for estimating fall harvest rates of 
females would be inappropriate because of the substantial 
natural mortality that occurs in the 7-9 months following 
tagging in late winter. 

We develop a joint known-fate tag-recovery model 
that uses radio-telemetry data to estimate tagging-to-har- 
vest survival and tag-recovery data to estimate harvest 
and annual survival rates. We use computer simulation 
to investigate the accuracy and precision of harvest rate 
estimates using this joint estimator as well as how the 
allocation of radio transmitters and tags affects preci- 
sion. We then illustrate use of this joint model for two 
species with different life-history characteristics: (1) 
female wild turkey with low tagging-to-harvest survival 
and low harvest rates, and (2) white-tailed deer with 
high tagging-to-harvest survival and moderate harvest 
rates. Also, for the white-tailed deer example, we show 
how the joint model can test whether individuals fitted 
with a radio transmitter have different harvest rates 
either because of hunter selection or behavior of the 
animal. 

Methods 

Kaplan-Meier known-fate estimator 

Radio-telemetry data can be used to estimate survival and 
harvest rates with the staggered entry design of the Kap- 
lan-Meier known-fate estimator (Kaplan and Meier 1958; 
Pollock et al. 1989). Survival can be estimated for prede- 
fined intervals and modeled as a binomial random vari- 
able and estimated by maximizing the likelihood of the 
product of interval likelihoods: 

g=l \ s / 

where 5^ is survival for each interval g (g = 1,. . .,G), dg is 
the number of individuals that died in interval g, and rtg 
is the number of animals fitted with radio transmitters 
that were being monitored at the beginning of interval g. 
The maximum likelihood estimate of an interval survival 
rate is 
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The product of survival rates from consecutive intervals 
is the cumulative survival rate: 

The staggered entry Kaplan-Meier estimator allows 
individuals to enter and leave the population throughout 
the study and assumes the following: (1) the sample is 
representative of the population, (2) survival of individu- 
als is independent of tagging, (3) survival of each individ- 
ual is independent of other individuals, (4) newly tagged 
individuals have the same survivorship as previously 
tagged individuals, (5) censoring, or removing the indi- 
vidual from the population available to be sampled, is 
random and not related to the fate of the individual (Pol- 
lock et al. 1989). 

The cost and logistics to affix radio transmitters and 
monitor a large number of individuals can be prohibitive. 
For example, the approximate cost of a VHF transmitter 
for ungulates is $250. A technician can reasonably moni- 
tor 40-50 deer with transmitters on a weekly basis, and 
each full-time technician may cost $20,000/year (C. S. Ro- 
senberry, PGC, pers. comm.). Deer fitted with Global 
Positioning System (GPS) transmitters do not need to be 
monitored by a technician, but cost $2,000-3,000. How- 
ever, if the GPS transmitters store information, as 
opposed to transmitting data via satellite or cellular tele- 
phone networks (an additional cost), they must be moni- 
tored and retrieved. Consequently, cost and logistical 
constraints often limit sample sizes, which can make it 
difficult to detect biologically important differences 
between groups. Powell et al. (2000) suggest >25 individ- 
uals in each marking group are necessary to make such 
inferences from a combined mark-recapture and radio- 
telemetry model. For the survival rate analysis. Pollock 
et al. (1989) recommend a minimum of 40-50 individu- 
als should have radio transmitters at all times, and addi- 
tional transmitters should be deployed prior to periods of 
high expected mortality. Winterstein et al. (2001) simi- 
larly suggested that 50 individuals were necessary to 
obtain a 95% confidence interval ±20% about an annual 
survival rate of 0.50. 

In addition, there is concern that the presence of a 
radio transmitter may influence a hunter's decision to 
harvest an animal (FuUer 1990; Etter et al. 2002; Mayer 
et al. 2002; Jacques et al. 2011) because radio transmit- 
ters can be visually detected by humans. The estimator 
may be biased positively or negatively depending on 
whether hunters avoid or select animals with radio 
transmitters for harvest. For white-tailed deer hunters, 
harvest decisions may vary according to characteristics 
of the deer, such as antler size and sex, and hunter 
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characteristics, such as age, experience, and geographic 
region (Jacques et al. 2011; PGC, unpubl. data). A con- 
founding cause for biased harvest rates could be that 
fitting animals with radio transmitters influences their 
behavior and makes them more susceptible to harvest 
(Caswell et al. 2012). 

Brownie tag-recovery estimator 

Dead-recovery models are special cases of capture-recap- 
ture models in which an individual can only be recap- 
tured once as a dead individual (Nichols 1992). The 
Brownie tag-recovery model is a parameterization of the 
tag-recovery probability when recovery is a result of hun- 
ter harvest, as described by Brownie et al. (1985), and 
can be used to estimate harvest mortality and annual sur- 
vival. An individual, upon being tagged, can survive the 
year (S), be killed by a hunter (K), or die from natural 
causes (I - S — K). If it is killed by a hunter, the indi- 
vidual could be retrieved (c) or not retrieved (1 — c). If 
retrieved, the tag could be reported (?.) or not reported 
(1 — l). In many studies, the probabilities of being 
killed, recovered, and reported are not separately estima- 
ble. The probability of being killed, recovered, and 
reported {Kc?i) is termed the recovery rate (/), and the 
probability of being killed and recovered (Kc) is usually 
termed the harvest rate (H). Consequently, when /i = 1, 
such as when sufficiently large rewards are used to 
encourage reporting of harvests, f=H (Nichols et al. 
1991; Pollock et al. 2001). We assume throughout the 
rest of this paper that rewards are used to ensure 1=1 
and/= H. 

Maximum-likelihood estimators for survival and har- 
vest rate are based on the number of individuals tagged 
and the number of tagged individuals reported in subse- 
quent hunting seasons (Brownie et al. 1985). We provide 
the likelihood for the Brownie model in the next section 
where we incorporate the known-fate data in a joint 
model. Tag-recovery models assume the following: (1) the 
sample is representative of the population, (2) age and 
sex are correctly determined, (3) no tag loss occurs, (4) 
survival of individuals is independent of tagging, (5) the 
year of a tag-recovery is correctly recorded, (6) the fate of 
each tagged individual is independent of the fate of other 
tagged individuals, (7) the fate of a tagged individual can 
be represented as a multinomial random variable, (8) all 
individuals in a given sex-age class have the same annual 
survival and recovery rates, and (9) annual survival and 
recovery rates can vary by year, sex-age class, and by area 
(Brownie et al. 1985). 

The primary advantage of tag-recovery studies is that 
they may be less expensive than other monitoring meth- 
ods, especially when the cost of capture is low and a large 
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number of individuals can be captured. Likewise, tag- 
recovery data may be more cost-effective when the cost of 
using radio transmitters to monitor fates of animals is 
high. Using reward tags can increase the reporting rate 
and improve the precision of estimates (Nichols et al. 
1991; Pollock et al. 2001). In addition, tags may be useful 
when a hunter's decision to harvest an individual is 
dependent on it being marked or unmarked, and the tags 
are not visible to hunters prior to harvesting an individ- 
ual. 

One disadvantage to tag-recovery models is that 
according to the parameterization of the model, all of the 
individuals that are tagged are considered available to be 
harvested at the start of the first hunting season after tag- 
ging. If mortality occurs between tagging and the first 
hunting season, the harvest estimator is negatively biased. 
Consequently, studies using tag-recovery models have 
attempted to capture animals immediately prior to the 



We use the subscript G to distinguish this tagging-to- 
harvest survival rate from the annual survival rate, S;, esti- 
mated in the Brownie model. 

The standard Brownie tag-recovery likelihood estimates 
fi and Si under the assumption that no mortality occurs 
between tagging and the first hunting season. However, 
the additional parameter Sa must be estimated if mortal- 
ity occurs during the tagging-to-harvest interval. The 
parameters Sq and /) in the following Brownie tag-recov- 
ery likelihood are confounded and not separately estima- 
ble, and this multinomial likelihood assumes that 
individuals were only reported in years when individuals 
also were released: 



L{Sci,Hi,Si\uii. . Mij) = ]^P(Mii,. 



P(uii, . . ., Uij) = — 



- SG,H,"«(SG;SiH,+i (Sg,S;S,+i. . .S,-iH,)""{l - 



Yluij]{ai-Ri) 



a,-R, 



hunting season. Another disadvantage of tag-recovery 
models is that cause-specific mortality cannot be deter- 
mined, although this is irrelevant if the objective is simply 
to estimate harvest rates. 



Joint Kaplan-Meier i<nown-fate and Brownie 
tag-recovery model 

Incorporating the Kaplan-Meier known-fate estimator 
into the structure of a Brownie tag-recovery model relaxes 
some assumptions of the Brownie model and can 
improve precision of parameter estimates. Animals fitted 
with radio transmitters are used to account for mortality 
between tagging and the first hunting season, such that 
animals fitted with tags do not have to be captured 
immediately prior to the hunting season. We illustrate the 
joint model using a monthly Kaplan-Meier known-fate 
model and a one-age-class Brownie tag-recovery model 
for a three-year study, but any study duration for estimat- 
ing tagging-to-harvest survival can be used. We used the 
Kaplan-Meier estimator as described in the previous 
section to estimate the cumulative survival rate between 
tagging and harvest, Sa, from the monthly survival rates, 
Sp, for year i, as 



where, 

fl; = the number of individuals tagged and released at 

the start of year i {i = 1, . . ., /), 

Uij = the number of individuals tagged in year i that 

are reported in year j ij = i, ...,/), 

Sg; = n ^ tagging-to-harvest survival rate in 

year i {i = 1, . . ., /), 

Hi = the harvest rate in year i {i = 1, ...,/), 

Si = the annual survival rate in year i {i = I, . . ., ]), 

\jji = SaHi + SciSiHi+i + ... + SciSiSi+i. . .Sj-iH/, and 

Ri = Uii + M,,i+i + ... + Uij. 

However, the known-fate estimates of survival are inde- 
pendent of the tag-recovery data so the combined likeli- 
hood for the joint estimator is simply the product of the 
two likelihoods: 



L{Sa,Hi,S,\dg 



n 



P{uii, . 



Sg; - W Sgi- 



The monthly survival rates estimated with the known- 
fate data allow Sg; and /; to be separately estimable when 
using the joint Kaplan-Meier known-fate and Brownie 
tag-recovery estimator (Table 1). 

The joint model can be extended to include data from 
animals fitted with radio transmitters in the estimation of 
the harvest rate. If we specify the binomial likelihood of a 
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Table 1. A three-year, one-age-class joint known-fate tag-recovery 
matrix sinowing tlie expected number of recoveries in year j of individ- 
uals tagged in year / (E[/?,y]) based on the number of individuals 
tagged in year / (A/,), harvest rate (H,), annual survival rate (S,), and 
cumulative tagging-to-harvest survival rate (5g/). 



Year 


Number released 








1 


2 


3 


1 


W, 


/V,Sg,Hi 


A/,SgiS,H2 


A/,SgiS,S2H3 


2 


A/2 




N2SG2H2 


N2SG2S2H3 


3 


A/3 






A/3SG3H3 



deer fitted with a radio transmitter being harvested (fr) in 
year i as 

where h, is the number of animals harvested and OT; is the 
number of animals available to be harvested, the maxi- 
mum likelihood estimate of the harvest rate of animals 
fitted with radio transmitters, for year i, is 

h, 

Hri — . 

ffli 

Again, because the harvest rate estimates for animals 
fitted with radio transmitters are independent of tagged 
animals, the combined likelihood is simply the product of 



transmitters (0, 10, 25, 50, 75, 100, 150, and 200 radio 
transmitters) and reward tags (50, 100, 200, 300, 400, 500, 
and 600 tags). We simulated a scenario in which tagging- 
to-harvest mortality was not accounted for by pairing 
each reward-tag allocation with 0 radio transmitters and 
constraining the monthly tagging-to-harvest survival rate 
to 1.0. We identified two plausible life histories of game 
species with different tagging-to-harvest survival rates and 
harvest rates to assess how harvest and annual survival rates 
affected precision of parameter estimates. Species A had 
a cumulative tagging-to-harvest survival = 0.90, harvest 
rate = 0.60, and annual survival rate = 0.30 and species B 
had a cumulative tagging-to-harvest survival = 0.69, harvest 
rate = 0.10, and annual survival rate = 0.60. In our simula- 
tions, we estimated annual harvest and survival rates even 
though these parameters were defined as constant, and we 
did not evaluate models that included age structure, 
although different harvest and survival rates for each age 
class can be incorporated into the model. 

We evaluated absolute and relative bias, in which the 
latter was calculated as absolute bias divided by the true 
value of the parameter, of the standard tag-recovery esti- 
mator by averaging the harvest rates over all years across 
simulations with 0 radio transmitters. We evaluated bias 
of the joint estimator using averaged harvest rates over all 
years for simulations with 10, 25, 50, 75, 100, 150, and 
200 radio transmitters. 

We evaluated the precision of annual harvest rate esti- 
mates, Hi, using the coefficient of variation (CV) defined as 



L{SGi, Hi, Hri, Si\dgi, Hgi, hi, Uii, 



n 

i=l 



P{Uii, 



the individual likelihoods. The likelihood for the joint 
model can be written as 

This model can take advantage of additional data from 
animals fitted with radio transmitters to estimate harvest 
rates by setting H,- equal to H,.,. Additionally, one can 
evaluate the assumption that animals fitted with radio 
transmitters have the same probability of being harvested 
as animals with tags by comparing models with and with- 
out Hi = Hri- 

Computer simulations 

We evaluated the model for different combinations of 
radio transmitters and tags using Program SURVIV 
(White 1983). We performed 1,000 simulations for 
combinations of numbers of animals fitted with radio 




X 100%, 



and the coefficient of variation of the root mean squared 
error (CV(RMSE)), which incorporated both precision 
and bias, relative to the known value Hf. 

CV(RMSE) = ' ' ' . X 100%. 

Hi 

To assist with interpretation, we graphed CV for the 
harvest rate for each year and species. In addition, we 
compared CV(RMSE) for each combination of radio 
transmitters to the scenario in which no radio 
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transmitters were allocated, which would be equivalent to 
a Brownie tag-recovery estimator. 

Case study: wild turkey 

The first case study uses data from female wild turkey 
where mortality between the winter capture period and 
the autumn hunting season introduces substantial bias in 
harvest rate estimates if the joint estimator is not used. 

Hen wild turkey were captured in Pennsylvania, Janu- 
ary to March, 2010-2012, in which all birds were fitted 
with an aluminum leg band (1242FR8A rivet band; 
National Band and Tag, Newport, KY) and a subset of 
birds were fitted with a Platform Transmitter Terminal 
(PTT) transmitter (North Star Science and Technology, 
King George, VA). The PTT transmitters allowed moni- 
toring fates of hens prior to the hunting season (Octo- 
ber-November) and estimating survival rates between 
month of tagging and the hunting season. Leg bands 
were imprinted with a toll-free number to report the 
recovery of the band along with wording "$100 
REWARD." We assumed retention of leg bands and 
reporting rate was 100% (Diefenbach et al. 2000). 
Because of the January to March capture period, not all 
turkeys needed to survive the same number of months to 
the hunting season; individuals tagged in January need to 
survive 9 months, whereas individuals tagged in March 
only need to survive 7 months. The joint known-fate 
tag-recovery model can account for unequal tagging-to- 
harvest mortality given that not all individuals are tagged 
at the same time. 

We modeled tagging-to-harvest survival by month and 
age class (adult or juvenile) but pooled data across years 
because of limited data for the juvenile age class. We 
modeled tagging-to-harvest survival for each age class 
because juvenile hens have lower nest-initiation rates and 
thus should be at lower risk of mortality. However, we 
combined age classes to model a single harvest rate 
because by the fall hunting season all tagged birds should 
have the same risk of being harvested. We developed two 
joint models, one in which harvest and annual survival 
rates varied by year (no. parameters = 23), and a second 
in which harvest and survival rates were constant over 
time (no. parameters = 20). We developed a third model 
where the tagging-to-harvest survival rate was 1.0, which 
is equivalent to a Brownie tag-recovery model that 
assumes no tagging-to-harvest mortality, to illustrate the 
effect of assuming no tagging-to-harvest mortality on har- 
vest rate estimates. We used Program SURVIV (White 
1983) to maximize the likelihood, obtain parameter esti- 
mates and standard errors, and compare models using 
AICc. The wild turkey data and SURVIV code can be 
accessed from Dryad (doi:10.5061/dryad.26n2q). 
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Case study: white-tailed deer 

The second case study uses data from white-tailed deer 
where the survival rate between tagging and the hunting 
season is >0.90, and we demonstrate how the joint model 
can be used to assess whether the presence of radio trans- 
mitters influences decisions by hunters to harvest deer fit- 
ted with radio transmitters. 

Pennsylvania has both known-fate and tag-recovery data 
for adult (>l-year old) and juvenile (<l-year old) female 
white-tailed deer for 2009-2011 in Wildlife Management 
Units (WMU) 2D, 2G, 3C, and 4B. Deer were captured in 
winter (January-April) when using bait is more effective 
because of limited food availability, but not harvested until 
autumn (October-January). Reward tags were small ear 
tags designed to mimic fur coloration by placing the white 
stud (6383 blank white stud; Destron Fearing Du™Flex Ear 
Tags, National Band and Tag Company, Newport, KY) on 
the inner part of the ear, and the black button retainer 
(6350 blank black button; Destron Fearing Du™Flex Ear 
Tags, National Band and Tag Company, Newport, KY) on 
the outer part of the ear. This design made it difficult for a 
hunter to identify the deer as tagged prior to harvest. Each 
tag was labeled with a toll-free telephone number, a unique 
identification number, and wording indicating a $100 
reward (Buderman 2012). To ensure retention of at least 
one mark until harvest, deer were fitted with a reward tag 
in each ear. We fitted deer with a VHF radio transmitter 
(M2510 with PET mortality; Advanced Telemetry Systems, 
Isanti, MN), GPS transmitter (G2000 Remote Release GPS; 
Advanced Telemetry Systems, Isanti, MN), or GPS/GSM 
transmitter (GPS Plus with Global System for Mobile Com- 
munications, Vectronic Aerospace, Berlin, Germany). 
Transmitters were labeled with a toll-free telephone num- 
ber to report the harvest, but no reward. We monitored 
survival weekly, and upon the transmission of a mortality 
signal from a radio transmitter, we investigated the fate of 
the deer. If a deer with a transmitter was harvested by a 
hunter and reported via the toll-free number, we recorded 
date of the harvest mortality. 

We first analyzed the known-fate data to identify the 
best model of tagging-to-harvest survival to simplify the 
known-fate portion of the model. When estimating 
monthly tagging-to-harvest survival rates, we considered 
deer captured at <12 months old to remain juveniles until 
June, which is the median date of fawn births (PGC, 
unpubl. data). We pooled data across WMUs because few 
deer died each month. We developed an a priori candi- 
date model set of how harvest and annual survival rates 
varied among groups. We structured the model to 
account for deer tagged in January that must survive a 
longer period of time from tagging-to-harvest than deer 
tagged in April. We always estimated different harvest 
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rates among WMUs because harvest management deci- 
sions were WMU-specific. The candidate set of models 
evaluated whether harvest rates differed by age class, 
among years, and between ear-tagged and radio-coUared 
deer. In cases where harvest varied by age or year, annual 
survival rates also varied by age or year. No annual sur- 
vival rate is calculated for radio-collared deer, so we did 
not model a tagging effect for annual survival. We used 
Program SURVIV (White 1983) to maximize the likeli- 
hood, obtain parameter estimates and standard errors, 
and compare models using AIC^. The white-tailed deer 
data and SURVIV code can be accessed from Dryad 
(doi:10.5061/dryad.26n2q). 

Results 

Simulation results 

Absolute bias in the Brownie tag-recovery estimator was 
less for species A than species B, reflecting the lower har- 
vest rate (Table 2), but relative bias was three times 
greater for species B compared with species A, corre- 
sponding to the lower tagging-to-harvest survival rate. In 
contrast, the joint estimator was essentially unbiased 
using both metrics, absolute and relative bias, for species 
A and B. (Table 2). 

As harvest rate decreased across example species, CV 
across all sample size allocations increased (Fig. 1). For 
example, precision across all years for species A ranged 
from 4 to 18% with an average of 8.8% across all simula- 
tions. Precision for species B ranged from 15 to 57% with 
an average of 27.8%, which was approximately three 
times as high as species A. With only a few radio trans- 
mitters to estimate tagging-to-harvest survival, precision 
of the joint estimator was poor (Fig. 1). This was most 
likely caused by increased model complexity, because only 
five parameters would have been estimated in the stan- 
dard Brownie tag-recovery model, whereas 32 parameters 
were estimated in the joint model. 



For both species, a feasible number of radio transmit- 
ters and reward tags would result in a lower CV(RMSE) 
for the combined model compared with the standard 
Brownie tag-recovery model. For species A, the joint esti- 
mator outperformed the tag-recovery estimator using only 
50 reward tags (the minimum tested in this study) and 
when at least 20 radio transmitters were allocated (Fig. 1). 
Using 100 reward tags and at least 25 radio transmitters 
was sufficient for the joint estimator to achieve a lower 
CV(RMSE) than the tag-recovery estimator for species B 
(Fig. 1). The difference between CV(RMSE) for the two 
estimators across the different allocations of reward tags 
and radio transmitters ranged from —3.22 to 6.61% for 
species A and from —13.11 to 17.47% for species B, 
where negative numbers indicate a larger CV(RMSE) for 
the joint model compared with the Brownie tag-recovery 
estimator. 

Case study: wild Turkey 

We captured and banded 162-169 adults and 74-128 
juveniles and fitted 32-55 adults and 6-15 juveniles 
with radio transmitters per year (Table 3). Survival dur- 
ing the January-September tagging-to-harvest period 
was 0.452 (SE = 0.0519, 95% CI = 0.35-0.55) for aduks 
and 0.514 (SE = 0.0906, 95% CI = 0.34-0.68) for juve- 
niles. The model with harvest and survival rates varying 
by year had the lower AIC^ value (166.3) compared 
with a constant model (AIC^ = 169.6). Harvest rates 
ranged from 0.012-0.045 when we only used data from 
leg-banded birds and did not account for mortality that 
occurred between tagging and the hunting season 
(Table 4). Consequently, not accounting for the tag- 
ging-to-harvest survival rate (approximately 50%) would 
result in harvest rates that were half that of truth. 
However, the coefficient of variations (CV) for harvest 
rate estimates for the joint model was similar to 
the model assuming no tagging-to-harvest mortality 
(Table 4). 



Table 2. Summary statistics of absolute and relative bias of harvest rates calculated from a standard Brownie tag-recovery model (assuming 
5g/ = 1) and a joint known-fate tag-recovery model based on computer simulations for specified tagging-to-harvest survival (Sg/) and harvest (/-/,) 
rates. 



Absolute bias Relative bias 



Model 




Hi 


X 


Minimum 


Maximum 


SD 


X 


Minimum 


Maximum 


SD 


Species A 






















Brownie 


0.90 


0.600 


-0.062 


-0.065 


-0.060 


0.001 


-0.104 


-0.108 


-0.101 


0.002 


Joint 


0.90 


0.600 


0.002 


-0.004 


0.010 


0.003 


0.003 


-0.006 


0.017 


0.005 


Species B 






















Brownie 


0.69 


0.100 


-0.030 


-0.032 


-0.026 


0.001 


-0.303 


-0.320 


-0.263 


0.014 


Joint 


0.69 


0.100 


0.002 


-0.002 


0.013 


0.002 


0.016 


-0.017 


0.126 


0.024 
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Figure 1. Three-year average CV for harvest rates from a joint known-fate tag-recovery estimator for (A) species A (cumulative tagging-to- 
harvest survival rate = 0.90, harvest rate = 0.60, annual survival rate = 0.30) and (B) species B (cumulative tagging-to-harvest survival rate = 0.69, 
harvest rate = 0.10, annual survival rate = 0.60) across radio transmitter allocations from 10 to 200 and three-year average percentage point 
difference between CV(RMSE) of harvest rates of the Brownie tag-recovery estimator and the joint known-fate tag-recovery model for (C) species 
A (cumulative tagging-to-harvest survival rate = 0.90, harvest rate = 0.60, annual sun/ival rate = 0.30) and (D) species B (cumulative tagging-to- 
harvest survival rate = 0.69, harvest rate = 0.10, and annual survival rate = 0.60). Negative values indicate a larger CV(RMSE) for the joint 
estimator compared with the standard estimator. 



Case study: white-tailed deer 

Between 0 and 69 deer were available in each month for 
each sex, age, and WMU group and between 0 and three 
deer died during each monthly interval. Annual WMU- 
specific availability of radio-coUared females over the 
hunting season ranged from 10-55 for adults and 0-15 



for juveniles, with <7 adults and <3 juveniles harvested 
each year in each WMU. If no radio-tagged deer were 
available during a harvest interval, harvest rate was con- 
strained to equal 0. For the Brownie tag-recovery portion 
of the analysis, 39-82 adult and 7-50 juveniles were 
tagged each year in each WMU, with 1-9 adults and <9 
juveniles harvested each year in each WMU. 
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Table 3. Number of adult and juvenile hen wild turkeys captured and 
fitted with $100 reward leg bands or radio transmitters, Pennsylvania, 
2010-2012. 



Year 


Adults 




Juveniles 




Banded 


Radio transmitter 


Banded 


Radio transmitter 


2010 


162 


55 


74 


6 


2011 


167 


42 


128 


11 


2012 


169 


32 


109 


15 



The best monthly survival model, obtaining 28% of the 
support, included only an age effect, with adult monthly 
survival equal to 0.9905 (SE = 0.0014) and juvenile 
monthly survival equal to 0.9799 (SE = 0.0048) across 
individuals from all WMUs. Because little mortality 
occurred between tagging and harvest, harvest rates from 
a Brownie model that was equivalent to the best joint 
model but did not account for tagging-to-harvest mortal- 
ity were underestimated by only 0.005-0.009. Support 
among the candidate model set for females was split 
between a model with no age, tag type, or annual varia- 
tion (AICc weight = 0.50) and a model with different har- 
vest and survival rates by age class (AICc weight = 0.45; 
Table 5). No models that included an effect of tag type 
on harvest had an AIC^ weight >0.01, such that we were 
not able to detect any evidence of radio transmitters 
affecting the behavior of hunters to harvest a deer or 
deer acting in a way that makes them more likely to be 
harvested (Table 6). 

Discussion 

To obtain accurate harvest estimates from a Brownie tag- 
recovery model for study designs where tagging-to-harvest 
mortality occurs, auxiliary information regarding tagging- 
to-harvest survival must be available. Otherwise, tagging- 
to-harvest mortality results in underestimation of the 
harvest rate. When the tagging-to-harvest survival rate 
was 0.90, harvest rate was underestimated by about 10% 
of the true parameter value and when the cumulative 



tagging-to-harvest survival rate was approximately 0.70, 
harvest rate was underestimated by 30% of the true value 
of the estimate. Whether natural resource managers 
should be concerned about potential bias depends on 
how underestimating harvest rates wiU influence manage- 
ment decisions. However, when tagging-to-harvest mor- 
tality data are available, they can be used in the joint 
known-fate tag-recovery model to minimize bias with 
little loss of precision. 

In the computer simulations, harvest data from radio 
transmitters did not contribute to harvest rate estimation, 
which is why there was little benefit to deploying more 
than 20 radio transmitters. If harvest data from radio 
transmitters were used to inform harvest rates, precision 
would have improved as more transmitters were 
deployed. However, whether this would be cost-effective 
would depend on study objectives and the additional 
costs of transmitters and personnel. We did not compare 
the joint model to a model with just radio transmitters, 
because the assumption of no tag-bias could not be 
addressed. However, if transmitters are small and not visi- 
ble, and tagging-to-harvest mortality is minimal, then 
harvest can be estimated directly from a Kaplan-Meier 
estimator. However, this may not be a cost-effective strat- 
egy if there is a high degree of mortality in the tagging- 
to-harvest interval. For example, to have a given number 
of turkeys with transmitters alive at the beginning of the 
hunting season, twice as many would have to be captured 
during the trapping season. 

A sufficient level of precision (CV) for parameters used 
to make wildlife management decisions has been 
described as <12.8% (Robson and Regier 1964; Skalski 
and Millspaugh 2002). However, the required CV may be 
dependent on the harvest rate of the organism and the 
management goals. Given a life history similar to species 
A, one with high harvest rates and low annual survival 
rates, the desirable CV range of 10-15% was attainable by 
allocating 10-20 radio transmitters and >50 reward tags. 
Conversely, 47 radio transmitters would have to be 
deployed to achieve a CV of 12.5% using only a Kaplan- 
Meier estimator for harvest rates. As harvest rate declines 



Table 4. Estimated han/est rates (H) for hen wild turkeys captured as juveniles or adults during January to March in Pennsylvania, 2010-2012, 
for a joint known-fate tag-recovery model that incorporated mortality between tagging and the first hunting season compared with a model that 
assumed no tagging-to-harvest mortality. 



Year 


Incorporating tagging-to-harvest mortality 




Assuming no tagging-to-harvest mortality 




H 


Sf(H) 




95% CI 


H 


SE{H) 




95% CI 


2010 


0.054 


0.024 


41.8 


0.01-0.10 


0.028 


0.011 


40.2 


0.01-0.05 


2011 


0.088 


0.0251 


28.5 


0.05-0.15 


0.045 


0.0116 


25.8 


0.02-0.07 


2012 


0.022 


0.0109 


49.5 


0.01-0.05 


0.012 


0.0059 


49.2 


0.00-0.02 



'CV = Sf(H)/H X 100%. 
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Table 5. Model selection statistics for joint known-fate and tag- 
recovery models of white-tailed deer harvest and annual survival rates 
for wildlife management units 2D, 2G, 3C, and 4D, Pennsylvania, 
2009-2011. All models separately estimated harvest and survival rates 
by wildlife management unit. The null model estimated harvest and 
annual survival rates by wildlife management unit but with no varia- 
tion among age class, year, or tag type. 



Variables 


Log-likelihood 


AAlCc 


AIC^ weight 


Null 10 


-328.7 


0.0 


0.50 


Age class 18 


-320.8 


0.2 


0.45 


Year 22 


-319.5 


5.8 


0.03 


Tag type 14 


-328.2 


7.1 


0.01 


Age and tag type 26 


-317.7 


10.2 


<0.01 


Age and year 42 


-305.8 


18.7 


<0.01 


Tag type and year 34 


-316.4 


23.7 


<0.01 


Age, tag type, year 64 


-294.7 


41.2 


<0.01 


'Number of parameters. 








Table 6. Estimated han/est rates (H) for ear- 


■tagged and radio-tagged 


female white-tailed deer for the best joint known-fate tag-recovery 


model that included a tag effect, Pennsylvania, 2009-201 


1. 


Ear tagged 




Radio tagged 


WMU H 


SE(f) 


H 


SE{f) 


2D 0.138 


0.0234 


0.145 


0.0447 


2G 0.116 


0.0248 


0.109 


0.0265 


3C 0.111 


0.0196 


0.150 


0.0461 


4B 0.133 


0.0209 


0.152 


0.0351 



to approximately 0.10, and tagging-to-harvest survival 
decreases to approximately 0.70, a CV of 10-15% for har- 
vest rate estimates was not achievable given less than 200 
transmitters and 600 reward tags. A number of possible 
sample size allocations, however, from 30 radio transmit- 
ters and 600 reward tags to 200 radio transmitters and 
350 reward tags could achieve a CV of 20%. In contrast, 
using only radio transmitters to estimate harvest rate 
would require 350 radio transmitters to achieve a CV of 
20%. We found, across both model species, only about 20 
radio transmitters were required for the joint model to 
achieve levels of precision equivalent to the Brownie tag- 
recovery model, which assumes no tagging-to-harvest 
mortality occurred. 

In cases in which the joint model is being employed to 
reduce estimator bias, CV(RMSE) may be a more represen- 
tative metric than CV for determining when using a joint 
model would be most beneficial. This was particularly true 
with the harvest rate estimates for species B, which would 
have exhibited more bias from tagging-to-harvest mortality 
than species A. By allocating at least 20 radio transmitters 
and 100 reward tags the joint model could improve CV 
(RMSE) by 0-15 percentage points for a species with low 



tagging-to-harvest survival (<0.70). In scenarios with high 
tagging-to-harvest survival, the benefit of incorporating 
radio transmitters to reduce bias may be negated by the loss 
of precision. For high harvest rates, 0.60 in the case of 
species A, the CV(RMSE) of the joint estimator was never 
more than twice that of the Brownie tag-recovery estima- 
tor. However, when the tagging-to-harvest survival rate was 
low, the reduction in bias was much more likely to warrant 
incorporating radio transmitters into the study design. 
Generally, only 20-25 radio transmitters were needed to 
improve CV(RMSE) of the joint estimator over the 
Brownie tag-recovery estimator regardless of the tagging- 
to-harvest survival of the study species. 

More complex models may be used to estimate harvest 
rates with the joint model than the simple model we evalu- 
ated via computer simulation (annual variation with no 
age or sex structure), but we would expect precision to 
decrease as model complexity increases. For example, 
including age structure would double the number of 
parameters (see white-tailed deer case study, Table 5). 
Another type of complexity that can be modeled with the 
joint known-fate tag-recovery model is differences in 
harvest rates for animals fitted with radio transmitters 
versus those fitted with simple tags. Most studies of harvest 
rates on wildlife species have monitored individuals with 
radio transmitters and assume that they are harvested at 
the same rate as individuals without transmitters. However, 
there is concern that the presence of a visible transmitter 
may affect a hunter's decision to harvest an individual 
(Fuller 1990; Etter et al. 2002; Mayer et al. 2002; Jacques 
et al. 2011) or that the transmitter leads to behavior by a 
marked individual that increases the probability it will be 
harvested (Caswell et al. 2012), leading to harvest estimates 
that are not representative of the population. Moreover, it 
is not clear how selective harvest by hunters might influ- 
ence harvest estimates because hunters could avoid harvest- 
ing animals with visible radio transmitters in the belief they 
are illegal to harvest or that harvest would adversely affect 
the research project (F. E. Buderman, pers. obs.). Con- 
versely, animals fitted with a radio transmitter could be 
viewed as a novelty and thus more likely to be harvested. 

To date, research to address the potential bias intro- 
duced by using radio transmitters to estimate harvest 
rates have used hunter questionnaires or simulated hunt- 
ing situations to assess hunter harvest behavior (e.g., Jac- 
ques et al. 2011). Such studies may provide insight into 
hunter behavior but are fraught with the difficulty of 
assessing whether a hunter's response to a hypothetical 
situation would be the same as their response to a real 
harvest opportunity. The model we developed is an 
opportunity to objectively assess the effect of transmitters 
on harvest by using animals marked with visible radio 
transmitters and less visible tags to separately estimate 
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harvest rates and statistically evaluate if parameter esti- 
mates differ. However, this method does not allow differ- 
entiation between hunters selecting for or against 
transmitters and marked individuals behaving in a way 
that makes them more vulnerable to harvest. If it were 
determined that there was no effect of the visible radio 
transmitter on harvest rates, as we found with female 
white-tailed deer in Pennsylvania, animals with radio 
transmitters can be used to increase the sample size (pre- 
cision) for harvest rate estimates. 
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